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IS.  ABSTRACT 


The  Fast  Fourier  Transform  (FFT)  and  modular  arithmetic  are  two  distinct 
techniques  which  recently  have  been  employed  to  increase  the  efficiency  of 
numerous  algorithm  in  the  area  of  symbolic  and  algebraic  manipulation. 

Motivated  by  work  done  on  fast  large  integer  multiplication  by  Schonhage  and 
Str  issen  [ LI]  and  by  Knuth  [7],  this  paper  analyzes  the  question  of  when  these 
two  techniques  can  be  utilized  concurrently.  The  desirability  of  the  convolution 
property  of  the  FFT  suggests  a  practical  definition  for  the  support  of  an  FFT, 
while  a  generalization  of  the  modular  rings  of  integers  motivates  a  resonable 
definition  of  a  finite  computation  structure,  A  Finite  Computation  Structure 
is  defined  to  be  a  commutative  ring  with  unity,  and  of  finite,  non- zero  charater- 
istic.  This  report  first  completely  characterizes  the  modular  rings  of  integers 
which  support  the  FFT  by  considering  the  prime  factorization  of  the  modulus. 

This  reoult:  Theorem:  Let  R  be  a  finite  computation  structure  of  characterization 
m.  Then  R  will  support  a  K-point  FFT  if  K  divides  p-l  for  each  prime  p 
dividing  m.  The  paper  then  concludt  s  with  examples  of  the  application  of  this 
result  to  the  problems  of  computing  products  and  powers  of  symbolic  multivariate 
polynomials. 
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Section  !_•  Introduction 

Ihe  Past  Fourier  Transform  (FFT)  ai  l  modular  arithmetic  are 
two  distinct  techniques  which  recently  have  been  employee  to 
increase  the  efficiency  of  numerous  a,  f>>rithms  in  the  erea  of 
symbolic  and  algebraic  manipulation.  motivated  by  work,  done  on 
fast  large  integer  multiplication  by  Set  nhage  and  S^assen 

and  by  Xnuth  [7],  this  paper  analyzes  tne  ^tlon  ™hilitv  of 
two  techniques  can  be  utilized  concurred  jy.  The  desirability  oi 
the  convolution  property  of  the  Fru  suggests  a  practical 
definition  for  the  support  of  an  FFT,  v  slle  a  generalization  of 

the  nodular  rings  of  integers  motivate!-,  a^rfS^^Comritition 
of  a  finite  computation  structure.  A  Finite  Computation 

Structure  is  defined  to  be  a  commutative  r Lng  with  unity,  and  oi 
finite,  non-zero  characteristic.  This  report  first  completely 
.  characterizes  the  modular  rings  of  intege  rs  which  support  the  t FI 
by  considering  the  prime  factorization  ox  the  moduius.  lhis 
characterization  is  then  extended  to  provide  the  following 
resSt:  Theorem:  Let  R  be  a  finite  computation  structure  of 

characteristic-'  m.  Then  R  will  support  a  K-point  PFT  if  *  divides 
p-1  for  each  prime  p  dividing  m.  The  paper  then  conclude.* 
examples  of  the  application  of  this  result  to  the 
computing  products  and  powers  of  symbolic  multivariate 

polynomials. 

Section  2.  The  Discrete  Fourier  Transform 


Definition:  Let  R  be  a  commutative  ring  with  unity,  written  as 
ifnr^Jntep  >  1,  and  WK  an  element  of  R  of  order  K?  ige^a 
primitive  K-th  root  of  unity.  Then  the  DISCRITE  PO.UFlMj, 
TRANSFORM  (DFT)  of  the  K-sequence  (a0,&1,..  .,a(K-i;j  is  tne  k 
sequence 

given  by  the  following  equations: 


i  =  0 


0<j<K-1  (1) 
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definition:  Assuming  the  same  conditions  as  above,  and  also, 
that  ^possesses  a  multiplicative  inverse  in  R  (i.e.  1/K)  then 
the  IWVggE  DISCRETE  Am  TRAKSTORM  (IDET)  of  the  K-kquence 
tau,a 1 , • . . ,a(E— l ))  is  the  K— sequence 


given  by  the  following  equations: 


°<i<K-1  (2) 


If  we  consider  the  term  ai*WK**(-j*i)  to  be 
equivalent  term 

•> 

v  «  ai*WK**((K-i)*j) 

then  we  can  rewrite  the  above  equation  as 


rewritten  as  the 


a. 

K-i 


i 

W 

K 


j 


f 


if  we  now  define  a(K)  =  a(Q),  then  the  inverse  DiT  can  be 
computed  from  the  DFT  by  merely  "flipping"*  the  input  sequence. 
?urf’  consists  m  replacing  the  i-th  term  by  the  (K-i)- 

th  term.  Thus  the  same  computational  algorithm  (for  the  DPT)  can 
be  used  to  compute  both  the  DPT  and  the  IDFT.  As  might  "be 
expected  from  the  terminology,  under  the  right  conditions,  the 
two  transforms  are  inverses  of  each  other,  and  thus  provide 
different  representations  of  K-sequences.  The  remainder  of  this 
F’pcr  ' wi. 11  be  concerned  with  determining  some  of  the  "right" 

?£ndl£iSn/  wh™?«ahe  DP?  Cfm  inverted.  Note,  also,  that 

£5?®  DPT  (and  the  IDFT)  are  linear  transformations  from  R**K  to 
R**K,  (where  R**K  is  the  ring  of  K-tuples  of  elements  of  R  with 
addlJ?on  and  multiplication)  since  the  quantities 
WK**(i*o)  are  all  "constants"  for  each  application  of  the  DITT. 
for  more  iniormation  on  this  approach  to  the  phenomenon  of  the 
DFT,  see  Nicholson  [10]. 

pie  computation  of  the  DFT  by  classical  techniques  usually 
involved  0(K**2)  operations,  as  the  computation  of  each 


_ ^ 
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transformed  element  took  K  multiplications  followed  by  K-1 

additions.  However,  Cooley  and  Tukey  [5]  demonstrated  a 
computation  scheme  by  ■  which  the  DPT  of  a  K-sequence  could  be 
computed  in  only  0(K  log  X)  operations.  This  method  has  become 
known  as  the  Fast  Fourier  Transform  (FIT).  The  key  concept  in 
the  reduction  of  the  computation  time  involves  the  fectofization 
of  K.  In  other  words,  for  K  highly  composite,  the  DFllcan  be 
computed  by  forming  sub-sequences  of  the  original  sequence, 

performing  a  DFT  of  fever  elements  on  them,  then  assembling  the 

resulting  sequences.  The  goal  here  is  not  to  develop  the  theory 

surrounding  the  FFT  (see  Cooley-Tukey  [5])  but  rather  to  instill 
some  feeling  for  the  immense  efficiency  of  this  algorithm  and 
thus  to  motivate  the  desire  to  find  as  many  ways  as  possible  in 
which  it  can  be  invoked.  During  the  remainder  of  this  paj^er,  the 
terms  DFT  and  FIT  will  be  used  interchangeably,  as  they  refer  to 
a  computation  function  and  a  computation  algorithm  which 
correspond  to  each  other. 

The  particular  virtue  of  the  DFT  in  many  applications  results 
from  the  followings 

Definitions  Let  A  =  (a0,a1,..  .,a(K-1))  and  B  =  (b0,b1,..  .,b(K- 
VT  two  K-sequences  in  R.  Then  the  CONVOLUTION  OF  A  AND  B, 
written  A*B,  is  the  K-sequence  C  «=  (cO,c1 , V.  .  ,c(K-1  )7~ where  tHe 
cj  ere  defined  as  follows: 


0<$CK-1  (3) 


where  bn  for  n  <0  is  defined  as  b(n+K). 

Convolution  Property  of  the  DFT: 

let  A  and  B  be  as  in  the  precedi..g  definition.  Let  A#.  and  B'  be 
the  DFT's  of  A  and  B  respectively.  Then  the  following  equation 
is  trues 

(A*B)'  =  A'  x  B'  (4) 

where  "x"  means  component- wise  multiplication  of  K-sequences.  In 
other  yords,  the  DFT  transforms  the  convolution  operation  in  R**K 
into  the  component-® ultiplication  in  R**K,  according  to  the 
following  commutative  diagram: 


PAGE 

t  - 


D FT  x  DFT 

R**K  x  R**K  — - - - >  *****  x  R**K 

a 


N*N  | MX" 


QED. 


As  a  result  of  (4),  the  convolution  of  two  K-sequences  can  be 
computed  in  the  following  way: 

1)  compute  the  transforms  A'  and  B'.of  A  and  B 

respectively. 

2)  perform  componentwise  multiplication  on  A  and  a  - 

to  obtain  a  K-sequence  C'. 

3)  perform  the  inverse  DFT  on  C'  to  obtain  the 

convolution  sequence  C, 

Thus,  the  DFT  provides  a  method  (though  not  an  intuitive  one)  for 
computing  the  convolution  of  sequences,  assuming  that  the  inverse 
DFT  of  step  3  is  possible;  i.e.,  if  we  can  turn  the  bottom  arrow 
around  in  the  above  figure.  What  this  requirement  amounts  to  is 
that  we  must  able  to  compute  a  "true"  inverse  DFT. 

Let  us  now  look  into  the  requirements  for  the  invertibility  of 
the  DFT.  If  we  apply  the  inverse  DFT  of  equation  2  to  the 
sequence  of  elements  computed  via  the  DFT  of  equation  1 ,  the 
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following  equations  .will  hold: 


■Jl) 


a  the  IDFT,  as  defined  originally,  will  be  a  true 

inverse  <*=>  a** hat-hat  sub  1"  =  a"sub  1"  for  all  1  <  K-*1  < _ > 


*9*  jU-JU  r 

4eO 


X  ik«i 


K«l. 

'  %  T<  toS 

x  77*0  < 


(Recall  that  delta(l.J)  *  1  if  iej  and  0  otherwise.)  Dotation: 
Inasmuch  as  the  following  quantity  appears  quite  frequently  in 
the  rest  of  this  papery  we  use  the  following  'notation: 

.X 

c  t  X  i  K-l 

Using  this  notation,  we  obtain  the  requirement:  The  IDFT  is  an 
inverse  <k=sk>  S“suh  1M  *  K  delta(0,l  ).  Thus,  a  reasonable  basis 

a  useful  DFT  (i.e.,  one  able  to  produce  the 
convolution;  is  to  define: 

Mln}tlop;  Let  K  be  >  1,  The  Support  of  a  K-point  DFT  is  n 

commutative  ring  R  with  identity  suHH%Jt““  ~  ~  ^ - 

(SI )  There  exists  ,1/K  in  R. 

;S?;  There  exists  wk  in  R  whose  oixier  =  K. 

K  delta(0,l)  for  ell  1  <  K-1. 

Note  that  the  support  of  an  FFT  is  an  “essential*  property  of  a 
ring  ;  i.e.,  if  a  ring  R  supports  a  K-point  ITT  and  R  is 

contained  in  a  larger  ring  R#,  then  R'  also  supports  K-point 


EFT's. 

Two  examples  of  support  of  a  K-point  DFT  are  t, 
with  wK  being  the  complex  number  cos(2WK)+ 
also  the  algebraic  number  field  Qfrj,  where  t 
solution  to  the  equation  x**K  =  1,  and'Q  is  the  ffi 
numbers. 


he  complex  number?; 
v- ■  sin(2%'K),  and 

is  a  formal 
eld  of  rational. 


Section  3«  Finite  Computation  Structures 


This  section  deals  with  the  problem  of  developing  a  sufficiently 
rich  generalization  of  the  structure  of  modular  rings  of 
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ttodSiie  to5dOTeloJVdatensSSct!li|EewE!^iz5ti?nJderives  froB 
is  precisely  modular^rithmetic  U^lrh£ri?iUllde£iyln£:  arithmetic 
Of  a  data  struck  “Jittin^*  3 15<£!?t"S?aliJ?s 

thinking,  we  propose  the  followings  BT  A"  a  result  of  this 

./IMte  Computation  Structure  is  a  commutative  ring 

th?  SlIite?i^f7^ffkS^f!?5!^e;1Bnc-  <Bc£S  that 
—  0  fnr  nil  a  -in  ^  18  ttk.  smallest  integer  m  such  xt — ■*- 

u  for  all  a  m  the  ring.  See  Lang  [6]  and  Albert  [1]?) 


the  _ 

ma  =  0  for  all 


that 


computation  structuret(fos)hofdcWacte^Lticam  It  iS  a  finite 

fsiSri^i  -“**ss»°s^rs[.a 

sssssr the  ciaas  ^ui”"Kof  ssg s?  concept 

integers  a. 

3)  for  811  n  811(1  all  n. 

4  Z  m,  ^rta  il?i^cSloi:  WhCTe  Euitiplic8tion  is 

p  uhere  Matipiication  ia 

(p7iifiMe£2  alrsf-5)  for  811  FrJ-s 


Section  4.  Characterization  of  Modular  Rings  of  Integers 


structures  tvis^vistthe2|prthit  CisSfirst  npp£ite  ^“putation 
the  case  of  aodSEr  riSL  ne£eSGary  to  consider 

rings  is  em^^in  e^°^nteg?Sf  one+?f 

StoT'-m26  those  modular  rings  of  infers  whiih  f  K- 

Klllrt  pi].1  ^his  vork1E'pi!l^  S?"  uE°nt®  down  by 
conditions  for  the  suroort  of  FPr'c;  ,fortlV  n<rcessary 

of  ^Pollard  :s  reLlt  "I The^ef'  *1 )  “el  imi  *?  °he  ""vlrie 

Pollard  (Theorem  C)  ’  ell*lnate  a  case  untouched  by 

the  resultingIequivilence^(!aeoremID)?V<>  SFtiFrSgFSZ,**? 
this  section,  the  author  assvsnes;  the  LadLP  S  ln 

understanding  of  elementary  number  theorv  UTin  gLS  hffi? 
and  elementary  group  theoxy  (as  in  Birkhof f  and  Mac-L^^p ] ) C  3 } 


Notes  In  the  remainder  of 
notational  conveniences: 


the  paper  w.i  will  use  the  following 
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Printed 

pi 

ei 

pi**ei 

Zpi**ei 

Z(m) 

Z(pi**ei) 

tlm) 

vK 


Written 

e.-P‘ 

c*’  e  ' 

-P; n 
7i/m9: 

*/p 

4>cW 


leaning 

p  sub  i,  a  prime  number 
e  sub  1,  a  positive  integer 
p  sub  i  to  the  e  sub  i  power 
The  ring  of  integers  modulo  m 
The  ring  of  integers  modulo 
the  ei— th  power  of  the  prime  pi 
The  multiplicative  group  cf  Zr 
The  multiplicative  group  cf  Zpi**ei 
phi  of  in,  the  Euler  phi-function 
K-th  root  of  unity 


let  n  be  odd  and 


Theorem  A:  Pollard's  theorem* 

r  I  ^ 

“of  ?totdi?fe 

support'll  K  f°r  Sl'{  then  the 


we  go  back  to  the 


|||>of:  In  order  to  prove  this  result, 
retirements  for  K-point  ITT  support.  ’ 

(of)  Eince  K  divides  each  pi-f  for  all  i,  then  hand  pi  are 
relatiifely  prcae  for  all  i.  This  implies  that  K  is 
relatively  priae  to  the  powers  of  the  pi's  and  alrn  to 
the  product  of  the  powere  of  the  pi's.  HeSSe  l  ?s 

z£r1Be  t0  ■*  W’iS  lmPlieE  that  K  possesses  ^ 

(Eg)  By  hypo^esis,  there  is  an  eleaent  wK  of  order  K  for  each 
(S3)  Ihecase  for  1^0  f^is  eBen,t  aiso  has  order  K  within  Za. 

hSrtly  £nl  1  trlvlal-  ^  Proof  for  1*0  depends 


Assume  the  above  hypothesis 
noint^^FT11  IK=^>**e  ^cr  some  prime  p. 


with  the  additional 
Then  Zm  supports  a  K- 


Proof  of  Lemma: 

*?£***’  K  and  p  are  relatively  prime. 
t-  rpsuit»  K  has  an  inverse  m  Zp**e. 

^  By  Zp**^theS1S’  th6re  exists  an  element  wK  of  order 

(S3)  Again,  the  case  1-0  is  trivial,  so  we  assume  1  is 

rero-  ,  Consider  now  the  expressions  for  S"sub  1"  wP 

Sde  SattheK^^r"SUb^“  (VKil>-  Clearly,  the  lift  haKd 
nfd  „  it?  is  zero  as  wK  is  a  K-th  root 

2**?*  tyu^  W  need  onlv  show  that  SMsub  1»  =  0  (modulo 

d°  Showing  that  wK-1^  (module  p) 

xcr  then  p**e  would  divide  S"sub  1,!.  Now,,  if  wK  =  i 


As  a 

K  in 

non~ 
V/e 


f  £  “eLc1  *3£S  *T.“)Snf 

5"  (^jpjssr  Mev:M  sg.s«r*si 
g  tSgsrRf  - 

»ssj* 

(modSo  piie)!  Ihus*  wK  *  1  (E0dul0  P)  ^  S"sut  1«  =  6 

QED  for  Lemma. 

nrnrt/SfiJS?  *w}thN  ^  proof  of  the  Theorem,  we  now  let  m  - 
prod ( pi  ei,i,1,r)  where  the  pi  are  distinct  crimes  rv  Tommc  a 

we  know  that  Zpi**ei  will  sutJDorf,  k—Ty^i  nt  ^  lemma  A, 

particular,  S-sSb  1"  =  0  (aS^ei)1?^  Sh  i  Bv  2S*  l\S 
Chinese  Remainder  Algorithm  (see  Lipson  m  foj  A“5l 
discussion  of  the  various  Chinese  Remainder1  "Alpo^ithSs)  on'11  th^ 

S(p^iPi1ierWUlJ  pi#*el?  S  that  S«sS  =*6  (SodSo 

proa^pi  ei,i,i  ,r; ; ,  i.e.,  (modulo  mi.  Thus  the  *Vrf*»o 

conditions  are  satisfied  and  Zm  supports  K-point  IPT's. 

QED. 

Kie  next  Theorem  presents  the  converse  to  Theorem  A  and  as  snoh 
11?  “?  C°mpJetely  the  odd  "“bers  which  caS  “pj“t  K- 


srts  K~point  FFT's  and 


m  = 


^eorem  B:  If  Zm  suppoj 

I>rod(pi**ei,i,  J,r),  m  odd,  then  K  divides  pi-1  for  »n  i  onH 
Zpi^eiX'  il”i!  Bent  VK  ln  2,0  SUCh  11,04  «K  iE  °f  ordeJ  h“n 

Proof:  Since  Zm  supports  K-point  FFT's.  then  hv  «;i  v 
Inverse  in  Zm  and  hence  K  and  m  are  relatively  crime  *  t+  *vvn 

a-v-i-'M-AstS*. !  » 

.  V* 

V  :Zn - . - >  TTZpi**ei 

*  4«l 

given  by  Wx)  =  (x  mod  p1**e1,x  mod  o2**e2  »  mrw*  »?>*•««.  n  * 

!(-nSi%®a5ir  SSt  JS&18.S 

^  ~~  &&+**"> 

Kerens  the  order  of  Z(pi*#ei)  is  given  by  *a|  ' 


I 


We  know  that  if  wK  has  order  K  in  Z(m),  then  K  must  divide  the 
order  of  Z(m).  Since  X  does  not  divide  any  pi,  then  K  divides 
the  product  prod ( pi-1  ,i,1,r).  However,  we  now  show  that  K 

divides  each  factor  pi-1*  By  the  3rd  requirement  for  K-pcint  FFT 
support  (S3),  we  are  guaranteed  that  S"sub  1"  =  0  (modulo  m  = 
pi**ei),  for  all  1  <  K-1  and  not  zero.  Since  the  pi**ei  are  all 
relatively  prime,  Then  S'*sub  1"  =  0  (modulo  pi**ei)  for  all  i. 
Now  consider  the  order  of  wK  in  Z(pi**ei),  say  K1.  Clearly,  K1 
divides  K.  If  K1  does  not  equal  K  then  SMsub  Kl"  ~  1+wK**K1+ 
\&**2K1+...*=1  +1+1... =K  (modulo  pi**ei),  since  the  order  cf  wF  is 
Kl.  But  K  and  pi*+ei  are  relatively  prime  for  all  pi  and  hence  K 
^  0  (modulo  pi**ei).  We  are  then  j  eft  with  the  conclusion  that 
if  Kl  ^  K,  then  S"a*sb  Kl"  ^  0  (mod  pi**ei),  which  is  a 
contradiction.  Thus,  wK  has  order  K  in  each  Z(pi**ei).  Thus,  K 

divides  the  order  of  the  group,  i.e.,  (pi-1)pi**(ei-l5.  Since  K 

and  pi  are  relatively  prime,  then  K  divides  pi-1  for  every  i. 

•QED. 

Having  finished  categorizing  completely  the  case  of  odd  moduli, 
we  consider  now  the  possiblity  of  even  mDduli  supporting  K-point 
ITT's.  This  possibility  is  discarded  by  the  following  theorem. 

Theorem  C :  No  even  modulus  can  support  an  EFT. 

Proof:  We  assume  here  that  K  is  non- trivial;  i.e.,  K  >  1*  In 

order  to  support  a  K-point  FFT  with  a  modulus  m  which  is  even,  we 
must  have  that  K  possesses  an  inverse  in  Zm.  This  condition  is 
equivalent  to  saying  that  K  and  m  must  be  relatively  prime  and 
hence  that  K  must  be  odd.  If  there  exists  an  element  wK  in  Zm  of 
order  K,  where  we  now  write  m  =  prod(pi+*ei,i?1,f)  with  pi  =2, 
then  by  means  of  the  isomorphism  described  in  Theorem  B,  the 
order  of  wK  in  Z(2**e1)  must  divide  K.  However,  since  the  order 
of  Z(2**e1)  is  2**(e1-1;.  the  order  of  wK  in  Z(2**e1)  can  only  be 
even  or-1.  It  can  not  oe  even  as  it  must  divide  K.  Thus  we  are 
left  with  the  fact  that  the  order  of  wK  in  Z(2**e1 )  must  be  1; 
I.e.,  wK  ss  1  (modulo  2**e1).  In  order  to  complete  our  proof  we 
look  at  the  possible  cases  for  m: 

Case  1;  m=  2**e1. 

In  this  case,  wK  has  order  1  in  Zm  =  Z2**e1  and  hence 
cannot  have  order  K. 

Case  2:  m  =  ml  *  2**e1  where  ml  is  odd. 

If  we  now  invoke  condition  S3  for  EFT  support,  then  S"sub 
1"  must  «  0  (modulo  ml  *  2**e1 ) .  Since  ml  and  2**d  are 
relatively  prime,  this  implies  that  S"sub  1"  =  0  fv/dulo 
2**e1)  for  all  1  <  K.  If  we  now  look  at  S"sub  *i"  v,wever, 
we  find  that  S"sub  1"  =  1  +  vK  +  vK**2  +  ...  =  -  +  1  +  1 

+  ..  ..  (modulo  2**e1)  (since  wk  =  1  (modulo  2*#e1))  =  K 
(modulo  2**e1),  which  is  not  equal  to  0,  as  K  is  odd. 
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ThUS’  S"EUb'1"  ls  not  zero  W»  have  a  contradiction. 

CTen°ir^dSi':^mot1lupp?rt0m^1Ctl0r'’  **  have  the  rpsolt  that 

Q5D. 

s«s istfszfi  r  r*«uits  of  this 

rings  of  integers  which  8u^?“SSt|Fi?aracterizatlon  of  modular 
i^r  liatf  pSTpfd^  Bf»'«  if  and  only  if  K  divides 

an.  eleaont  vK  in  Zm  of  o?der  K  fi 
The  argument  proceeds  as  follows  i 


2. 


Evei?  possess  primitive  roots  of  unitv 

whicA  LZthlUUiZ 

The  element  wi  defined  by 

£S»  W^SKSSS,  fiji . 

«  -  C12#-) 

3-  Wf&SSS  f«ny  of  Congruences"  (“  »> 

wK  =  wi  (modulo  pi**ei)  i  =:1,...,r 

»d'ChSeofnortcrTfnr  I?  Z(pl**el)  for  H  * 

element  of  theCesired  or^heeljLhflds’  S'?  caI1'>J.vBys  find  an 
existence  of  K-th  roots  ofSitv  iS  1th®  hyp°thesis  of 

Thoreirj  A,  B  and  C  is  redundSt.  1  ^  equivalence  Implied  by 

primes^rom^eoreSrA^d^  I^mSh  £  Jh®  Jro?thesis  Tor  odd 
aach  i  and  that  K  is  ntflf*  SS"g  S&Ert  SSV£l 


PAG!  12 


’ equivalence  of'ra???1?  °£  r^riotions  as 
exactly  the  statement  of  Theorem  D„  Thcore,flS  A»  B  and  C  result  In 


QED. 


q  v  * 

seef^n,'  wi^have^pfovided  presented  in  this 

nn^s  of  integers  supporting  th^Past^^i^r  ltodular 

P?Ytin»W1  +***  • 


Examples; 

SimiSrlvr^i.fnly  2'1?int  OT's  (with  w>:  =  -1) 

O-point  Jft's  eanCbe  performed°in  «?>  for  all  e. 


k  • 
2. 
3. 


oint  TTT's  can  be  r»rfn™£*  1  r^irit  "T's  for  a] 
Of  the  form  prSdf^ftSJ^2?  f?r  “1  * 
e  ivi-Vhe?S  +  '  itX’1  arbitrary 
e,g.,  17**©,  33**e',  17**e  *  33*ie% 


ggetion  5=  ^tension  to  Unite  Cogitation  Structures 


As  we  have  seen  in  np^tinn  a 

bbppprt  K-point  JIT's'.  le  !ov«ffi  5S^ai"  modular  rings  will 
of  the  finite  computation  structures.  io  result  t0  include  - 


some 


Theorem ;  If  r  is  . 

characteristic  m,  then  JR  sun™??  computation  structure  of 
for  every  prime  p  div?dIn|  ^ItE  a  K-point  HT  if  K  divide!  p!i 


5®0r?n  J  we  know^hat  lymZ^pSts5!1?6™?,  “•  Then  by 

the  requirements  SI,  S2  £nd  S3Pof a  result; 
5  Pcs  of  characterise?*  +?Ctl8n  2  hola  111  Z/mZ.  Sincp 

R  isomorphic  to  Z/mZ.  ^Becaus^Sl  thS?®  ?Jists  a  subring  S  of 
support  of  an  FFT  involve  onlv  r^°f  the  requirements  for 
isomorphic  copy  of  Z/m 2  will  satLf??:h??e^ati?ns,  then  any 
Thus  S  supports  JT-point  IW'*1  wL?he  5uPP°rt  requirement, 
support  is  essential^d^th^/^sTp^la  £i nlfe"  * 


Note;  If  Q-;Z/mZ  _ >  r  sc  • 

then  0*(wK)  is  a  K-th  root  of  ijn?fi0rp?isS  onto  a  subring  of  R, 
modular  K-th  roots  of  unity  will  be  LfttcLt^omPH^! 
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Section  6:  Applications 


The  results  presented  have  been  applied  to  the 
computing  products  and  powers  of  symbolic 
polynomials  over  the  integers  (see  [3]  and  [4]). 
brief  example  examines  the  application  of  the  FFT 
computation  structure  to  the  computation  of 
polynomj  als. 

Let  f  and  g  be  univariate  polynomials  over  some  coefficient  ring 
R,  with  the  degree  of  f  =  m  and  the  degree  .of  g  =  n.  We  can  thus 
write: 


problems  of 
multivariate 
The  following 
in  a  finite 
products  of 


+  ***** 


■£1*3=  K  t  +  •••  t  1 L*" 

^  -  b  o  +  b,  y-  +  bt  tl  ♦  •  •  b  u  x H , 

The  product  polynomial  h(x)  =  f(x)*g(x)  can  then  be  written  as 

h  (.*>  m.  t0  C(*  4»  Ct4t*  + 

where  the  coefficients  ,ci  are  computed  by 

P  •  ZL  /  ,  A  ,  U  , 

*  I.  D  ' 

If  we  set  K  ~  B+n+1  and  form  £  K-sequences 

A  =  l*o,  o)  ■ 

B  -  <•  b.  ,  !»,,  •  •  •,  VjM|  O,  ■ 

then  the  K-sequence  C  given  by 


0} 


C,  z.  (i,  t,  c.  ^ 

is  precisely  the  convolution  of  A  and  B,  as  defined  in  Section  2. 
As  a  result  of  this  observation,  if  the  coefficient  ring  R 
supports  K-point  FFT's,  the  the  product  polynomial,  h(x)  can  be 
computed  using  PFT's  as  follows: 

1)  Form  the  2  K-sequences  A  and  B 

2j  Apply  a  K-point  FFT  to  each  sequence,  yielding  A  and  B 

3)  Multiply  A  and  B  component-  wise,  yielding  in  C. 

4)  Apply  the  inverse  K-point  .FFT  to  obtain  C,  which  is 

the  sequence  of  coefficients  of  the  product 
polynomial  h. 

In  particular,  if  the  coefficient  ring  is  the  integer,  we  cannot 
perform  FIT  s  as  there  are  no  roots  of  unity  other  than  1  and  -1. 
However,  if  the  largest  coefficient  (in  absolute  value)  in  A,  B 
or  C  is  bounded  by  M,  then  we  can  embed  f(x)  and  g(x)  in  the 
polynomial  ring  Zm[x],  where  m  >  2M  and  Zm  supports  K-point 
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fcIpSte  tte’1productPPllCatl0n  °f  the  4  steps  Itentioned  atK  ve  “H 
h  11O  E  >0 

However,  since  thn  maximum  coefficient  is  less  then  m/2,  then 
h(x)  =  f(x)g(x)  ana  we  are  finished.  '  ’  en 

Similarly,  if  the  coefficient  ring  R  =  Z[y]  (i.e.;  f  arx  g  are 

f3Vanc&tr  lrVX  y.,over  the  integers),  we  center  bed 

i  and  g  in  “[£»y]  where  m  is  sufficiently  large  and  Zrfvl 

SS2rtfi?”p?1?i  *** f  •  ^in»  application  of  the4  steps  g?ien 
above  will  yield  a  polynomial  1  1 

gc-u,  ^  c  vy 

f(x.y)fi(x?y)!  16  Ch0sen  aPPr®Pri-c«iy,  h(x,y)  actually  equals 

pie  computation  of  powers  of  polynomials  proceeds  in  an  aralopous 
manner.  otudies  performed  on  these  two  computational  rroblems 
have  indicated  that  the  ITT  methods  provide  significant 
improvement  m  the  efficiency  of  these  algorithms  fofaSe 
class  of  polynomials.  (See  Bofmeau  [3]  [4]).  E 


Section  7s  Summary 

presented  several  results  relating  modular 
.?che“efi  811(1  the  last  Tourier  Transform.  In 

Sr1  JnSJSSi66?  0f  rln£s  of  Regers  in  which  the 

of  computed  is  completely  characterized  by  the  prime 
decomposition  of  the  modulus.  Also,  an  extension  of  this  result 
■^or  computation  structures  similar  to  modular  rings  of  integers 
yields  a  sufficiency  hypothesis  for  the  computatiSf  of  m”  ^ 

Further  research  in  this  area  might  be  motivated  by  the  desire  to 

fin?+«neCeESa+y+^ 811(1  sufficient  conditions  for  TFT  support  in 
finite  computation  structures  or  in  general  rings.  1  As  an 
Nicholson  [10],  using  the  definition  of  FFT  to  be  a 
transformation  which  takes  "convolution-product*1  rings  into 

su?,S?^t‘"fr02Vcr  •  ringf ’  Provides  the  result  that  m* s  are 
if f ln.-Ttfion  liags  (i.e.,  rings  with  no  zero  divisors) 
11  nJ-r  v  rln£  contains  a  particular  primitive  root  of  uni  tv.  It 

ri^eSlrable  t0  extend  “is  type  of  ^ult  to  sparger 
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